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1.   Introduction  and  Summary 

Estimates  of  location  parameters  based  on  rank  order  tests 
In  the  univariate  one-sample  and  two-sample  problems  have  been 
considered  by  various  authors.   For  the  details  of  references, 
the  reader  Is  referred  to  the  most  recent  papers  of  Gastwirth 
[4]  and  Bloch  and  Gastwirth  [J>].      In  this  paper  we  consider  the 
estimates  of  location  parameters  In  the  multivariate  one-sample 
and  two-sample  problems.   These  estimates  are  obtained  through 
the  use  of  coordinate  wise  rank  order  test  statistics  such  as 
the  Wllcoxon  or  normal  scores  statistic  considered  by  the  authors 
In  [9]  and  [11].   The  distribution  of  these  estimates  is  shown 
to  be  symmetric  with  respect  to  the  parameters  being  estimated. 
These  estimates  are  also  translation-invariant  and  robust. 
Furthermore  the  estimates  are  shown  to  be  asymptotically  normal. 
Finally  the  asymptotic  efficiencies  of  these  estimates  relative 
to  the  estimates  based  on  the  vectors  of  means  and  medians  are 
obtained  by  applying  the  criterion  of  Wilks  generalized  variance 
(cf.  Anderson  [1],  p.  166). 


2.   Point  Estimates 

L^t  ^a  =  t^la--"V^'  a  =  l,...,m  and  Yp  =  (Y^^, .  .  .  ,Ypp) , 
P  =  l,...,n  be  independent  vector  values  random  variables 
drawn  from  p-variate  absolutely  continuous  cdfs  (cumulative 
distribution  functions  P(x)  and  G(x)  =  P(x-A)  respectively 
where  x  =  (x^,...,x  )  and  A  =  {A,,...,A  ).   For  each  coordinate 
consider  the  statistic 

where  E^^^  is  the  expected  value  of  the  a-th  order  statistic 

of  a  sample  of  size  N  =  m+n  from  a  distribution  ^.;  zi^'  =  1 

if  the  a-th  smallest  observation  from  the  combined  sample 

{X.T,...,X.  ;  Y. ,,..., Y.)  is  an  X  observation  and  Z,t   =  0 
Jl      jm'   Jl'     jn'  Na 

otherwise,  a  =  1,...,N;  j  =  l,...,p.   Let  for  each  j  =  l,...,p 

(2.2)  A*  .  sup{Aj  :  hjtXj^,...,X.^,  Y.^-A.,...,Y.^-Aj)  >;lj} 

(2.3)  Aj»  =  inf  {Aj  =  h3(Xj,,...,Xj^;  Yj^-Aj, . . . ,Yj„-Aj )  <  v,] 


where  p..  is  a  point  of  symmetry  of  the  distribution  of  h.  when 
J  J 

A.  =  0,  j  =  l,...,p.   It  was  shown  in  [6]  that  the  estimate 

J 

A„  .  =  (A*  +  At*)/2  of  A.  has  more  robust  efficiency  than  the 

n  m 


classical  estimate  >    Y  .a/"^  -   >    X.  /m  in  the  corresponding 

univariate  problem.   Thus,  as  a  natural  extension  of  [6],  we 
propose 


.A 


(2.^)  A„  =  {A„,  ,  .  .  .  ,A,,  ) 

^    '  ~N    ^  Nl'    '  Np' 


as  an  estimate  of  A  for  suitable  functions 


(2.5)  h  =  {h^,...,h    )    . 


Now  let  us  consider  the  one-sample  problem.   Suppose 
Z  =  (2,  , ...,Z   ),  a  =  1,...,N  is  an  independent  sample  from 
a  p-variate  cdf  F(x-  9)  =  F{x, -0,  , . . .  ,x  -9    )  where  P  is 
absolutely  continuous  and  is  diagonally  symmetric  about  zero. 
(Let  us  recall  that  P  is  diagonally  symmetric  about  zero,  if 
its  density  remains  invariant  under  simultaneous  changes  of 
signs  of  all  the  coordinate  variates  (cf.  [11]).)  As  before 
consider  for  every  univariate  sample  X^*^^  =  (Z  .-,,...  ,Z  .„) , 
a  statistic 

(2.6)  hj*  =  hJ(z'J')  =l!^E^J'  Z^'/N,     J=1,...,P 

where  E*,^'  is  the  expected  value  of  the  a-th  order  statistic 
of  a  sample  of  size  N  from  a  distribution  ^.(x)  given  by 

J 

(2.7)  ^*(x)  =  ^**(x)  -  Y**(-x)   for  x>0  and  ^*(x)  =  0,  otherwise, 

J        J        J  J 

Y.  (x)  is  symmetric  about  zero  or  uniform  over  (-1,1,),  and 

J 

Z^^^   =  1  if  Z .^  >  0  and  Z^^^  =  0  otherwise;  a  =  1,...,N; 
j  =  1, . . . ,p.   Denote 


(2.7)     9j  =  sup  le.  :   h*tzj,-ej,...,zj^-ej)  ^,.^5 


(2.8)     «r  =  inf  fSj  :  h;(Zj,-e.,...,Zj„-ej)  <^jj 


Where  ix  .  Is  a  point  of  symmetry  of  h*  when  9     =  0.   Then  for 
J  J       J 

suitable  functions  h.,  we  propose 


A  A  A  A  ^        ^^ 

(2.9)  ^N=    ^%1"'-'W'      ^NJ  =  ^^j  +^j  ^/2  ^    j  =  l,...,P 


as  an  estimate  of  9. 

/v  A 


The  estimates  A„  as  well  as  9     defined  above  form  a  general 
class  of  estimates.   Two  Important  members  of  this  class  (to 
which  we  shall  pay  specific  attention  in  Section  5  when  we 
discuss  the  relative  asymptotic  performances  of  the  estimates) 
are  (1)  the  Wllcoxon  type  estimates  ^^tj^\   and  ^^rj^\    resulting 
from  the  Wllcoxon  statistic  by  taking  for  f .  and  ft  respectively 

J  J 

the  rectangular  distribution  over  (0,1)  and  (11)  the  normal 

A  A 

scores  type  estimates  A„,-xx  and  £jt/Xa  resulting  from  the  normal 
scores  statistic  by  taking  for  T.  or  '^**   the  standard  normal 
distribution.   The  present  authors  investigated  the  properties 
of  the  tests  based  on  the  statistics  (2.1)  and  (2.6)  for  the 
p-variate  one-sample  and  c-sample  problems  in  [11]  and  [9] 
respectively.   This  paper  examines  the  properties  of  the  estimates 
A^,  and  9^,  based  on  the  statistics  (2.1)  and  {2.6)    respectively. 


5.   Regularity  Properties 

The  following  theorems  are  Immediate  extensions  of  the 
theorems  proved  by  Hodges  and  Lehmann  [6]  for  the  case  p  =  1 
and  are  therefore  stated  without  proofs. 

THEOREM  5.1.   The  distribution  of  the  estimate  A„  as  well 
as  that  of  9 „  Is  (absolutely)  continuous  If  P  Is  (absolutely) 
continuous. 

THEOREM  3.2.   Invariance  Properties. 

(3.1)    A,T(Xn  ,.  .  .,X  ;  Y.+a,.  .  .  ,Y  +a)  =  A.JX,  ,  .  .  .  ,X  ;  ¥,,...,¥  )  +  a 
^-^   ■'    /vN^  1'    '  m'   1  ~*    '  n  ~'    ~N^  1'    '  m'   1'    '  n'   ~ 

{5-2)   £n^^i+^-"'V^'  =  In^^i-'-'^n^ +^ 


where  a  =  ( a^ , . . . , a  )  Is  any  1 x  p  vector  of  constants. 

THEOREM  3.3.   The  distribution  of  the  estimate  A„  defined 


In  (2.4)  Is  diagonally  symmetric  about  A  If  either  of  the 
following  two  conditions  hold; 

(a)  The  distribution  F  Is  diagonally  symmetric,  and  for 
each  j  =  1 , . . .  ,  p 

(3.3)   h.(X.,  ,  .  .  .,X.  ;  Y.T,...,Y.  )+  h.(-X.^  ,.  .  .,-X.  ; 

j^  jl'      jm'   jl'      jn'    j'   jl       jm 


-Y.-,  ,  .  .  .  ,-Y.  )  =  2|i  . 

jl'    '   jn^    ^j 


and 


(3.^)  h.{X.^+Si,...,X.^+a;   Y  .-^+a,  .  .  .  ,Y  .^+a) 

=  hj(x.^,...,Xj^;  Yji,...,Yj.^+a  for  all 


(b)  The  sample  sizes  m  and  n  are  equal  and  for  each 
j  =  1 , . . . , p ;  h  .  satisfies 

(3.5)   h,(X.,  ,.  .  .  ,X.  ;  Y., ,...,¥.  )  +h.(Y.-,  ,  .  .  .,Y.  ;  X.-,,...,X.  )  =  2|i.  . 

and  {5^i<^). 

THEOREM  3'^'      The  distribution  of  9^^  defined  In  (2.9)  Is 
diagonally  symmetric  about  9   If 

(i)   F  Is  symmetric  about  zero 

(11)   hJ(Z^Jb  +  hJ(-Z^Jb  =  2Hj;  J  =  l,...,p 

Remark:   The  conditions  on  h.  and  h.  In  Theorems  J>.J>   and  J>.h 

J      J 

are  satisfied  both  for  the  Wllconxon  and  Normal  Scores  statistics. 
THEOREM  5.5.   For  every  vector  a  =  (a, , . . . ,a  ) 


«V 


(3*6)   P[h  .  (X  .,,...  ,X.  :  Y  .T -a  ., .  .  .  ,Y  .  -a.)  <  u, . ;  j  =  l,...,p] 

<  P[A-,  <  a]  <  P[h.(X.,  ,.  .  .  ,X.  ;  Y  ., -a  . , .  .  .  ,Y .  -a.)  <  \i.', 
_   L^N   A,-"  —   I-  j^  jl'    '  jm'   jl   j'    '  jn  j^  —  ^j' 

j  =  1 , .  .  . ,  p  ] 
and 

(3.7)   P[h!(Zj^-aj,...,Zjjj-aj)  <  Hj;  j  =  l,...,p]  <   P[9jj  <  a] 
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4.   Asymptotic  Normality 
Let  m(N)  and  n(N)  for  N  =  1,2,...  be  a  sequence  of  pairs 
of  sample  sizes  tending  to  infinity  in  such  a  way  that 
m(N)/N  — >  A,  say,  and  let  A^  be  a  sequence  of  the  values  of 
the  parameter  A.   Also  for  the  one-sample  problem  consider  the 
sequence  of  sample  sizes  N  =  1,2,...  and  let  9     be  a  sequence  of 
values  of  0.      In  both  the  cases  we  shall  indicate  the  dependence 
of  h.,  h.  and  p..  by  writing  h^.,  h^ .  and  [i^ .   respectively.   Let 
us  denote  the  marginal  distribution  of  Y.   (or  Z.  )  by  F.(x) 
and  that  of  iY  .^,Y^^)    (or  (Z.^,Zj^^))  by  F.^^{x,y)   when  A  =  0 
(or  0=0)  and  assume  that  the  marginals  are  absolutely 
continuous.   Next  let  us  define 


1  1  2 

J    Jj(x)dx   -    fj    Jj(x)dxj        if    j    =  k; 


(4.1)        d.^. 


0 
oo      oo 


-co    -00  b  0 

if    J   ^  k;    Jj.    =  Y"^,    j    =   l,...,p    . 


00 

(4.2)  ^j    "/     {dJj-[Fj(x)]/dx]   dFj(x)    ;      Jj    =  ^ 


-1 

J 


t^-"     '}.^ 


-00 

^    f    J*   (x)dx     if   j   =  k 
0 

00         00 


'^ 


iff    jJ[2Fj.(x)-l]J^[2Fj^(y)-l]dFj^j^(x,y)      if   j  j^  k 


-OD       -00 


J  J    =    ^J""^     ,     j     =     l,...,p 


CD 

(4.4)     g.    =fldJ*[2F.{x)-l]/dx]<3iF.{x)    ;  jj  =  ^J"^ 
0 

Finally,  let  a  =  (aT,...,a^)  be  a  ix  p  vector  of  constants  and 
'>-'     J-      p 

let  c-,  , . .  .  ,c^,  . . .  be  real  constants  and  let 
(^•5)  Aj^  =  -  a/cj^   or   9^  ^   -  a/c^ 

Then  to  prove  the  asymptotic  normality  of  A^  and  9^   we 
shall  make  use  of  the  following  two  theorems,  the  first  of 
which  Is  an  Immediate  consequence  of  Theorem  6.1  of  Purl  and 
Sen  (1966)  and  the  proof  of  the  second  Is  given  In  [11]. 

THEOREM  4.1.   If  h.,  j  =  l,...,p  are  given  by  (2.1)  with 
^  . ,  j  =  1 , . . . , p  satisfying  the  assumptions  of  Theorem  6 . 1  of 

Purl  and  Sen  (I966)  and  If  m(N)/N  — >-  A  as  N  ->.  od  ,  then  the 

1/2 
random  vector  N  ^  (h„.  -IJ^vr-j  J  =  !>•••»?)  has  asymptotically  a 

p-varlate  normal  distribution  with  mean  vector  zero  and  covarlance 

CD 

matrix  D  =  ((dj^^))  given  by  (4.1).   [Here  Hj^^.  =J     Jj[Fj(x)  ]dP^.  (x) .  ] 

-00 
THEOREM  4.2.   If  h*   j  =  1, . . . ,p  are  given  by  (2.6)  with 

J 

T*  satisfying  the  assumptions  of  Theorem  4.1  of  Sen  and  Purl 

1/2  *  * 

(1967),  then  the  random  vector  N  '  (^m-  "^tj*-'  ^   ~   !>•••■»?) 

has  asymptotically  a  p-varlate  normal  distribution  with  mean 
vector  zero  and  covarlance  matrix  F  =  ((f.^))  given  by  (4.3). 

OD 

^^NJ  =/  j;[2Fj(x)-l]dF.(x).] 
0 


Remark:   (l)  The  limiting  distribution  of  N"'"/^(h  .  -  n   , 
j  =  l,...,p)  is  nonsingular  if  and  only  if  the  functions  J.  and 
F.  are  such  that  the  moment  matrix  of  J.[F.(x)],  j  =  l,...,p 
is  nonsingular.   It  is  singular  if  and  only  if,  a.s.  F 

(4.6)       Jj[Fj.(x)]  =  211  "i^  F^(x)  +  constant  [cf.  [9]] 


kT^J 


(a,  are  some  constants) 


(2)   The  limiting  distribution  of  N-"-/ ^(h*  -  p.   ,  j  =  l,...,p) 

is  nonsingular  if  and  only  if  the  moment  matrix  of 

(J .[2F . (x)  -  1] ,  j  =  l,...,p}  is  nonsingular.   It  is  singular, 

if  and  only  if,  a.s.  F 


(4.7)   j!(2Fj(x)-1)  =  ^^  \  J^[2F^[y)  -  1]   +   constant  [cf.  [11] 

(a,  are  some  constants)  . 

The  following  theorems  give  the  relationship  between  the 
limiting  distributions  of  h,,  and  A,,;  and  h,,  and  9    . 

THEOREM  4.3.   Let  G  be  the  continuous  p-variate  distribution 
function  whose  marginals  have  locations  zero  and  scale  parameters  1 , 
and  suppose 


^""'^^       ^ii"^  W%--^Nj^  ^^j  '   J-i'---p] 


=  G 


Z^V^^A  ^p^  Vp- 

^1  ^p 
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where  P^  Indicates  that  the  probability  Is  computed  for  the 
parameter  values  Aj^  or  9^   and  where  h^^ .  stands  for 
^NJ^^J1---'V(N^  V^"'V(N))  andh^.  for  h;;.(Z. ,,..., Z.^). 
Then  for  any  fixed  A  and  9, 

(4.9)  llm  PaIVAn  -  ^^  -^}  =  ^(-XT'  ••-  -H) 
N — >.  00  ~  ^1  P 

(4.10)  llm  ^el^An-  V  -^]=''(^"--'\^ 

N-— ?.oo  -v,   ''  ''  VI  p 

Proof:   By  Theorems  3.2  and  3 A 


llm  P.  fC„(A,T  -  A)  <  a\ 
N — >-  00  -v 

=  llm  P  S  C„  A-  <  a  ] 
..   ^  oC  N-vN  —  ^J 


a .  a . 

.  y i       Y   ,  , i)  <  LL 


^'"oo^a^^J^^Jl""'^J^^N)'  ^Jl-C^'  ••"^0n(N)-C^^  -Hj    ' 

J  =  1 , .  .  .  ,  pj 


N-> 


=   llm  Po{hj^^jr---'^jm(N)5  Yjl'""Yjm{N)5-^Nj  -°  ' 
N — >■  00  ~ 

J  =  1 ,  .  .  .  ,  p  j 

/  a,  Bn       a  B  •■ 
_  p  /  1  1        P  P 

\   ^1  ^P   . 

This  completes  the  proof  for  (4.5),  and  that  for  (4.6)  Is 
completely  analogous.   Thus  combining  the  results  of  Theorems 
4.1,  4.2  and  4.3  we  obtain 
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THEOREM  4.4.   (a)  Under  the  assumptions  of  Theorem  4.1, 
N  '    (A„  -  A)  has  asymptotically  a  p-variate  normal  distribution 
with  mean  zero  and  covarlance  matrix  x  =  (t.,  )  given  by 

d 
(^.11)  -^jk  -   X(li)e.e,  >         jMc  =  l,...,p 

1/2  '^ 
(b)  Under  the  assumptions  of  Theorem  4.2,  N  '  (0»t  -  ^)  has 

asymptotically  a  p-variate  normal  distribution  with  mean  zero 

and  covarlance  matrix  /\    =    (A.,  )  given  by 

(^.12)  A^j^  =  fj^i^/gjgi^  ;  J.k  =  l,...,p 


Special  Cases.   A.  Let  J.  be  the  Inverse  of  the  standard  normal 

J 

distribution,  and  J*  be  the  inverse  of  the  chl-distributlon 

J 

A  A 

with  one  degree  of  freedom.   Then  the  estimates  A„  and  9     reduce 
to  the  normal  scores  type  estimates  A„/Tr%  and  ^vrrX'v  respectively. 

In  this  case  the  covarlance  matrices  x  =  (x.,  1  and  A  =  (A.,  ) 

I        ,1   ,  .1        ,1  , 

reduce  to  x  =  (t.,)  and /Y  =  l^-k'  respectively  where 

(4.13)  T   =  Ji^-= —  ;     j,k  =  l,...,p 

J^   A(1-A)B  (F,i)B  (F,f) 


and 


I  .     ,$ 


(^.1^)  ^jk  =  ^^^-  ^^"^jk  5         j,k  =  l,...,p 

where 
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(^.15)    7ji,(F,l)  = 


if  j  =  k 


00   oo 


-OO  -00 


if  J  7^  k 


(4.16) 


B  (F,;^)  =  f    — — 4 


f .(x)dx 


-00 


(|)[I-^(Fj(x))] 


J  -  1, 


where  <|){')  is  the  density  of  the  standard  cumulative  normal 
distribution  function  ^(  •  )  •   B.  Let  J.  or  J^  be  the  Inverse  of 
the  rectangular  distribution  over  (0,1).   Then  the  estimates 
A„  and  9„  reduce  to  the  Wllcoxon  type  estimates  A„/t3\  and  d..fD\ 
respectively.   In  this  case  the  covarlance  matrices  t  =  ("^n-v) 


R 


R 


and  A  =  ^^-k^  reduce  to  x^  =  (t"j^)  and  A^  =  (^^1^)  respectively 
where 


(4.17)  T 


R 
jk 


[12A(1-A)B^(F,R)]-^ 


If  j  =  k 


00    00 


/   /  Fj(x)Pk(y)^j,kf-'y)-¥ 


^  -00   -00 


-1 


[A{1-A)B  {F,R)Bj^(F,R)] 


and 

(4.18) 

^jk  =  ^ti-^'-5K 

where 

CD 

B^.(F,R)    =    f     f|(x)dx    ; 

(4.19) 

-00 

Now  denote 

j  ,  k  =  1 ,  .  .  . ,  p 


J  =  1 , .  .  .  ,  p 
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where 

(4.21)  Aj^j  =  median  (Y^^,  .  .  .  ,Yj^)  -  median  (^ji^  •  •  •  »X  .j^^)  , 

J  =  1, • • • ,P 
and 

^Mj  "  median  (Z^^, .  .  .  ,Zjjj)  ,    j  =  l,...,p 

f^-22)     Yj^-Xj^=  (Y^-X^,...,Yp-Xp)  ,    Zj^=  (Z^,...,?p) 

where 

—   —     n  m  _     N 

In  the  next  section  we  shall  investigate  the  performance  of  the 
estimates  A„  and  9^  with  respect  to  the  estimates  (4.20)  and 

(4.22)  based  on  the  vectors  of  medians  and  the  means  respectively. 
To  this  end  we  shall  need  the  following  well  known  result. 

THEOREM  4.5.   ( a )  Under  the  assumption  that  the  variance 

of  marginal  distribution  function  F .  exists  and  is  finite  for 

& J  

each  j  =  1 , . . . , p , 


(4.24)    lim  P.[N-'-/^(Y.,-  X„)  <  u]  -  ;^  {Ut,...,u  ) 


where  ^        is  the  cumulative  distribution  function  of  the 
p-variate  normal  distribution  with  mean  vector  0  and  covariance 


Ik 


matrix  y~  =    ((cr.j^))  and    C.^  =   cov  (Xj^,Xj^^)  when  A  =  0  or  0  =  0. 

{i|.25)     llm  Pe[Nl/^(Z  -0)  lu]  =  5       (u-^,...,u  ) 
N-^  00  ~  lO,2_  J 


(b)  Suppose  the  marginal  distribution  F.(x)  Is  absolutely 
continuous  for  each  j  =  1 , . . . , p  and  Its  derivative  at  0  denoted 
by  f  .(0)  exists  and  Is  finite.   Then 

(4.26)    llm  P.[N^/^{L,- A)  <  u]  =  f  (u.  ,...,u  ) 


*  ** 


where  y~  =  ((  <^^-^))    Is  given  by 


I If  J  =  k 


(4.27)   a^jk=<^ 


l^otXja  "  0'  \p  ^01-^1  /fj(0)fi,(0)    If  J  ^  k 


and 


(4.28)     llm  P.[N^/^(0  -£)  <  u]  =  $     ^  (u   ...,u  )  . 
I^->GO  -^      ^  ^  -^      [0,7"]  ^  P 


5.   Asymptotic  Efficiency 

To  obtain  an  Idea  of  the  relative  performance  of  one 
estimator  with  respect  to  another,  we  employ  the  notion  of  "the 
generalized  variance"  a  concept  Introduced  by  Wllks  (see  Wllks 
(1962),  p.  547).   The  "generalized  variance"  of  a  p-varlate 
random  vector  (X-,  ,  .  .  .  ,X  )  with  non-singular  covariance  matrix 
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^  =  {p.^(y.<r^)    is  defined  to  be  var  X  =  (T-^    *  '  •  o'p  ^^^  "Pik" 
where  "det"  denotes  the  determinant. 

Suppose  that  the  two  asymptotically  unbiased  estimates 
T  and  T'  with  9   with  asymptotically  non-singular  matrices 

m        m     mm         rn  ?        ml    mi    m  t 

Y_     =   ^Pjk^'j  ^k^  ^"^^  ^   =  ^Pjk^j  ^k  ^  ^^^^i^e  N  and  N' 
observations  to  achieve  equal  asymptotic  generalized  variances. 

Then  the  asymptotic  relative  efficiency  of  T  with  respect  to  T' 

is  defined  as 

_-l 


(5.1)     e„  „,  =   lim 
^'^  N— >  00 


11 
N 


T' 


.T 


^i  •••  '^p 


det  II  p 


T 


T      T   ^  ^-  II  T  I 
0-1  •••  ^p   ^et  llPikl 


Now  from  (4.11),  (4.12),  (4.24),  (4.25),  (4.26)  and  (4.28)  the 
generalized  variances  of  A^^,  |^,  Y^-  X^^,  Zjj,  %^   and  0^  are  given 

by 


(5.2) 


d.  . 


var 


(N^/^A^)  .fy-^-Iii-.detllp^^i 


j=l  e"A{l-A) 

J 


where  d .,  and  e.  are  given  by  (4.1)  and  (4.2)  respectively,  and 


(5.3) 

(5.4) 
where 
(5.5) 


A     ,   /^l/2  .1/2 


j ,  k  =  1 , .  .  .  ,  p 


var  (N-^/^  e.,)  = 


l^^r../,f)ae^\\p%\\ 


J-k 


f .,/f 1/2  f 1/2 

Jk^  J  J   kk 


j,k  =  1, 


,P 


l6 


2 
(5.6)     var  [N^/^{Y^-X^)]  =  J^  ^^^  •  det  H  p  .  J| 


[A(l-A)]-P  var  (N^/2  Z^) 


(5.7) 


where 


(5.8) 


var  [N^/^  Ajyj] 


M 


j=l  i^A(l-A)f^{0)         J^ 
=  [A(l-A)]-P  var  (N^/^  9^) 


h[?AZ,^   >   0,    Z,_  >  0)  -i]    If  j  ^  K 


1    If  j  -  k 


ka 


Hence  denoting  e^^j,*  and  e^  ^*   as  the  efficiency  of  T  with 
respect  to  T*  for  the  one-sample  and  the  two-sample  problems 
respectively,  we  obtain 

-1 


(5.9)   e 


(2) 

(A„,Y„-X„) 
*  'V.N  'V.N  ~N 


I  (T.-  e  ,/d , , 


.J=l 


J   J   JJ 


det  lip   II  /det  llp^j^ll 


iP 


-1 


and  similar  type  of  expressions  for  e^^-'  _  ,  e^^''  ^ 


and  e^^^  /^ 


Our  main  interest  is  to  study  the  relative  asymptotic  performances 

A  A  f\  —     — 

of  A^w^r^  J  A^w„x,  \„   and  Y,,  -  X„  in  the  two-sample  case  and  those 
~N(^)   'vN(R)   -vM     <\.N  ~N 


of  9 


,    9     and  Z,,  in  the  corresponding  one-sample  case. 


^N,J'  -N,R'  -M  ^"""^  ^N 
[At  this  stage  we  would  also  like  to  draw  the  attention  of  the 

reader  to  a  paper  of  Bickel  (1964)  who  earlier  proposed  the  estimate 

A 

9-,T,_^  of  9   and  investigated  its  asymptotic  performance  with  respect 
'^N I K  J     'V' 

A  A 

to  the  estimates  9     and  Z„.   Since  the  estimate  0^T/^^^  is  only  a 

'uM     '\jN  •vNtRJ 

particular  member  of  a  relatively  more  general  class  of  our 
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estimates,  some  of  our  results  duplicate  his  results.   However 
this  duplication  is  retained  here  In  order  to  provide  compre- 
hensive comparison  and  for  the  ease  of  reference.] 
Thus  applying  the  definition  (5.I),  we  have 


(5.10)      e 


(2) 


=  e 


(1) 


^N{J)'~N"~N  ~N{^)'~N 


Act?  b?(f. 


•1    r 


J=l 


3    r-'^) 


-1 


det  lip 


Jk' 


_det  II  p 


i 


jk  ^ 


where  B.{P,J)  Is  given  by  (4.16)  and 


(5.11)   pi 


/^        00         00 


-00      -CD 


If  J   =  k 


(5-12)      e 


(2) 


-1^) 


~n(r)'Xn  ~N         ~N{R)'~N 
p 


r  120-^(Bf(P,R)) 


^p-^ 


det  ijpj^lj  /  det  jlp^j^ll  J 


-1 


where  B.(F,R)  Is  given  by  (4.19)  and 

J 


r  00   00 


12 


(5.13)  P^, 


/   /  Pj(^)I'k<y)<'Fj^^(x,y)4 


^  -00   -00 


If  j  7^  k 


If  j  =  k 
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(5.1^) 


.(2) 


-  Jl) 


^N®'%N(R)     §N®'N(R) 


T^   B^(F,$)/12B^(F,R) 


-1 


det  lip 


^  I 


det  II  p 


I 


jk  J 


(5.15) 


,(2) 


=  Jl) 


~N 


($)'^    ~N(J)^-M 


j=l  4f^(0) 


detllp^^ll 


det  II  p^  II 


-1 


where  p.,  is  given  by  (5.8). 


When  p  =  1,  the  efficiencies  {5-11),  (5.12),  (5-13)  and 
(5.14)  reduce  to  the  familiar  expressions  in  the  one-sample  as 
well  as  two-sample  univariate  case  [cf.  [5]^  [6]].   We  shall 
study  the  above  efficiencies  in  different  situations. 

Case  1.   Total  Symmetry.   Let  us  assume  that  the  estimators 

^N($)'  ^N(R)'  h'   Xn-In'  %,l'    %,R'    %  "-""^   ^N  ^^^^  Pairwise 
independent  coordinates.   In  this  situation  the  covariance 

matrices  reduce  to  diagonal  matrices  [cf.  [9],[11]]  and  we  have 


,(2)         _e(l) 

~n(r)'Xm"Xn   %(R)'^N 


(5.16)    e.^"'   _  _  =  e\^'  =   12 


fj.^. 


-1  -1 


U=i 


b^(f,r) 


(5.17)   e^     _  _  =  e. 


.?. 


'N(i)'^N  ~N    -N(^)'% 


2  ^2, 


TTaf  Bt(F,® 


j=l 


-1 
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(5.18)   el^)   ^ 


^N{J)'^N(R) 


e(2) 


1 


(5.19)   el^L  _  =  el^L 


=  4 


^'Xn~~N    ^*~ 


~N 


.J=l 


U=i 


J  '>^ 


T^  B^AF,f)Bf{F,R) 


-1 


-1 


Hence  from  the  results  of  Hodge s-Lehmann  (I961),  Chernoff -Savage 
(1958)  and  Mlkulskl  (I963),  we  have 


(5.20) 
(5.21) 
(5.22) 
(5.23) 


Inf  e^    _  =  0.864 
Inf  e^    _  =  1 


inf  e^ 


TT 


inf  e^  _  =0.35 


where  ji    Is  the  set  of  all  totally  symmetric  absolutely  continuous 
p-variate  distributions. 

Now  suppose  that  the  components  of  P  are  totally  symmetric 
as  well  as  identically  distributed,  then  the  efficiencies  are 
Independent  of  p,  and  hence  the  same  as  in  the  case  of  one- 
sample  univariate  situations. 

Case  2.   Equally  Correlated  Distributions.   Let  us  now 


assume 


that   the   distribution  of    {xy\xx^h    is   independent   of   J 


a      '    a 
and  k  when   0=0.      Then,    it   turns   out   that 


f\j  f\j 
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(5.24)      e^  _     =  e,         _    _ 

B  7  AY    -Y 

~N(R)''^N  ~N,R'~N  ~N 


=   12  Cr^  B^(F,R) 


l  +  (p-l)p 


12 


l  +  (p-l)p 


12. 


-1   r 


1  -p 


12 


1-  p 


12- 


1-p 


-1 


(5.25)      e. 


=  e. 


.? 


."^-..-x 


•N(5)'^N  ^N(^)'~N"^N 


of  B^(F,J) 


1+  (P-I)p 


-I    -1 


12 


1+(P-1)P*;, 


1-  p 


12 


1-p 


12-J 


1-p 


-1 


(5.26)     e. 


=  e. 


~N(J)'%,R  ^N($)'-^(R) 


,-2, 


Bj(F,f)B^''(F,R) 
1^ 


1+  (p-l)p 


*  7_-l   r 
12 


l  +  (p-l)p** 
12- 


1-  p 


12 


1  -p 


12-1 


1-p 


-1 


(5.2?)      e^     _     =  e^ 


=4^  fi(o) 


1-p 


12 


1  -Pi2. 


1-P 


-1  r 


1+  (p-l)p^2' 


r 


Ll+  (p-l)p'^'2- 


i-p 


-1 


Case   J>.      Normal   Case.      Let   the  underlying  distribution 
function  F  be  a  non-singular  p-variate   normal   distribution  function 


with  mean  vector  zero,    and   covariance  matrix  5~"  =    (cT.,  cr.o",   )    then, 
from   (5.12),    (5-10),    (5-14)    and    (5-15)    we   obtain 


(5.28)  e 


^N(R)^~N 


=  e 


-N(R)'~N"~N 


3 

IT 


det  II  p 


i_-l 


jk' 


det  Up 


Jk' 


where 
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00         00 


(5.29)  p;,=/      /    $.(i]  l^[il    df^^Jx.y) 


-00      -00 


(5-50) 


=   e. 


=   1 


'N(f)'XN  ~N         ~N,;|'^N 


(501)        e^ 


P.T.rrT\>§, 


N(J)'%{R) 


5 


det  \\p*^\ 


1    -1 
P 


det  |!pjj^„j 


=  e. 


A 


N($)'^N{R) 


where  p^,^  Is  defined  by  (5-29), 


(5.32)        ej^^        .     =  e.  =  | 

^N{R)'^  ~N{R)'-Sm 


det  IIp^JI 
detllp^f^l 


^P-^ 


Now  Blckel  (1964)  has  proved  that  for  p  >  2, 


(5.53) 


inf  e^    _  =  0 


^   -^N ,  R  ~N 

where  (j)  is  the  set  of  all  non-singular  p-variate  normal 
distributions  and,  since 


(5.34) 


inf  e. 


F^^  £n(5)'^n 


=  1   for  all  p 


it   follows   that 


(5.35) 


sup  e^  ^ 

F&'t>      £n($)''^N{R) 


=    OD 
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Thus  we  find  that  when  the  underlying  distribution  function 
Is  normal,  the  multivariate  normal  scores  estimator  A„/■x^  or 
0  .-r.  (which  is  as  good  as  the  means  estimator)  can  be  Infinitely 
better  than  the  multivariate  Wilcoxon  type  estimator  A„,J^^  or 

for  p  >  2.   This  leads  to  the  question  of  examining  the 


;ion 


~N(R) 

relative  performance  of  these  two  estimators  viz.   £j^CK)  i°^   ~N(¥)  ^ 

A  /\ 

and  0^t/t.\  or  (A,w„x)  for  the  case  when  the  underlying  dlstribut' 
function  F  is  blvariate  normal  NCOjCr^jCr^^'P)  •   The  efficiency 
behavior  of  O^j^j^   and  1^^^^^^  (or  A^^j^  and  1^^^^)    is  given  by  the 
following 


THEOREM  5.1.   The  efficiency  of  P^^pj  Mi  £n(?) 
and  A-jf-^i  is  Independent  of  cr,  and  a^   and  is  given  by 

-Tl/2 


^  ^N(R) 


(5.36) 


-N(R)'-N($) 


IT 


2 
2 


_i-_e: 


1-9(1-1  cos-1  |)' 


(1+2  cos  -I  (Tr+u))'Il/2 


u(Tr-u) 


where  u  is  determined  by  p  =  2  cos  (u+7r)/3-   The  function 
e^  is  monotone  increasing  for  0  <  p  "^  1  and  (2) 


A 

0 


N(R)'~N($) 
symmetric  about  p  =  0  and  hence  unimodal.   Finally. 


(5.37)    11m 


1/2 


1pI-^i  ^n(R)'%($)   ^ 


V   Sin  -^■ 
5 


=  0.91 
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The  proof  is  the  same  as  in  [2]  (see  Theorem  5.2  [2]) 
since  when  the  underlying  distribution  is  normal. 


©A     /^     =  e. 


~N(R)'~N{5)     ~N(R)'~N 


In  a  similar  way  the  efficiency  behavior  of  0„/■r•^  and  0,. 

A  A  "^ 

or  Ajj/^j   and  ^  may  be  obtained.      The  details  are  omitted. 
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